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The present work considers energy harvesting by implementing both thermo- and piezo-electric power 
generation modules on a bifurcating tube, which produces dual-temperature thermoacoustic oscilla¬ 
tions. The present system distinguished from the conventional standing-wave one does not involve heat 
exchangers and uses two different energy conversion processes to produce electricity. To measure and 
analyze the sound waves generated, an infrared thermal imaging camera, hot wire anemometry, and two 
arrays of K-type thermocouples and microphones are employed. It is found that the total electric power is 
approximately 5.71 mW, of which the piezo module produces about 0.21 mW. It is about 61% more than 
that generated by a similar conduction-driven thermo-acoustic-piezo harvester. In order to gain insight 
on the heat-driven acoustic oscillations and to simulate the experiment, thermodynamic laws are used to 
develop a nonlinear thermoacoustic model. Comparison is then made between the numerical and 
experimental results. Good agreement is obtained in terms of frequency and sound pressure level. 
Finally, Rayleigh index is examined to characterize the conversion between thermal and sound energy. In 
addition, energy redistribution between different thermoacoustic modes is estimated. It is found that 
lower frequency thermoacoustic oscillations are easier to trigger. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

Fossil fuel exhaustion and greenhouse gas emission are 
increasingly urgent issues. However, waste heat generated by 
burning fossil fuel is released to the atmosphere in industry, which 
is harmful to the environment and energy wasting. To recover 
waste heat, there is a need for energy conversion techniques. 
Thermoacoustics [1-6] is a research field which deals with the 
conversion from thermal to acoustical energy [7,8]. For refrigera¬ 
tion, heat should be transferred into the working gas of a ther¬ 
moacoustic system during the expansion [7], and dumped from the 
gas during compression. A thermoacoustic heat engine (also known 
as prime mover) operates on a thermodynamic cycle in the direc¬ 
tion opposite to that of the refrigeration cycle. Theories and models 
have been developed and reviewed for thermoacoustic phenome¬ 
non and devices [2,5,7]. 

Thermoacoustic engines are more attractive than the traditional 
Stirling engines [9 involving pistons, due to the fact that they lack 
moving parts, require little maintenance and last long. However, 
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these engines 10,11 provide lower efficiency, little reliability, and 
high fabrication cost. For example, the traveling-wave engines [11] 
have high efficiency and high reliability, but have high fabrication 
costs due to the presence of regenerator and heat exchangers. The 
standing-wave engines are reliable, and cheaper to fabricate; but 
they are inefficient. Typically, a thermoacoustic engine is one of the 
thermofluidic oscillators in which thermodynamic oscillations are 
generated and sustained by external temperature differences [12] 
achieved by heat exchangers that play critical roles in producing 
such oscillations. However, implementing such heat exchangers 
causes not only additional energy consumption but also increased 
manufacturing cost. 

Thermoacoustics systems involve the energy conversion from 
heat to sound or other form of mechanical work [13-15 . Optimi¬ 
zation of the energy conversion process was studied by several 
researchers [6,15 . Various forms of input energy have been tested 
such as solar power 13] and concentrated laser pair [16 . A mini¬ 
ature thermoacoustic engines were developed [17 . Recently, a 
thermo-acoustic-piezo system was designed and tested to achieve 
energy conversion from heat-to-sound-to-electricity [18 . It was 
shown that the output power was very low (approximately 
0.128 mW) and two heat exchangers were used. 
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In this work, a thermal energy harvesting system is designed and 
experimentally tested. For this, a cylindrical tube splitting into two 
bifurcating daughter branches is designed. As a heat source is placed 
inside the mother tube, it provides a mechanism to produce both ‘hot’ 
and ‘cold’ thermoacoustic oscillations (also known as Rijke-Zhao 
oscillations). To achieve energy harvesting in the Y-shaped tube, a 
piezoelectric diaphragm is implemented to the bifurcating branch 
with ‘cold’ oscillations. And a thermoelectric power generation 
module is attached to the end of the other bifurcating branch with 
‘hot’ oscillations. The experimental setup is described in Section 2. 
Thermodynamic measurements are then conducted. This involves 
measuring the temperature and monitoring the sound fields, as 
described in Sect. 3.1. The performances of the piezo- and thermo¬ 
electric modules are then measure and compared, as described in 
Section 3.2. Finally, in Section 4, thermodynamic laws are used to 
develop a nonlinear thermoacoustic model to gain insight on the 
heat-to-sound energy conversion and to simulate the experiment. 
Comparison is then made between the numerical and experimental 
results. In addition, Rayleigh index and sound energy redistribution 
are examined to characterize the energy conversion, and the energy 
redistributed between different eigenmodes. 

2. Experimental setup 

A non-axisymmetric bifurcating tube with a thermo- and piezo¬ 
electric diaphragm attached is designed and tested to demonstrate 
the possibility of harvesting thermal energy by using two different 
energy conversion processes, as show in Fig. 1. The angle between 
the bifurcating branches is 2a and 0° < a < 90°. As a heat source 
(flame or solar heater et al.) is placed inside the mother tube of 
length 0.3 m, it provides a mechanism to produce dual¬ 
temperature thermoacoustic oscillations [19,20 . It is worth 
noting that the oscillations generated in the present system are 
convection-driven ones, which are different from the conventional 
conduction-driven standing- or traveling-wave prime movers [1,7]. 
T-shaped junctions protruding from the bifurcating branches are 
designed for placing thermocouples, microphones or hot wire 
anemometry. Before the thermodynamic measurements are con¬ 
ducted, the microphones are calibrated by using a piston phone 
(124 dB at 250 Flz). The measurement error is about ±0.1 dB. The 
hot wire anemometry is calibrated by using a Pitot-static tube. The 
thermocouples are connected to a Pratos K-type digital 



Fig. 1 . Experimental setup of dual-temperature thermoacoustic system with a piezo¬ 
electric and thermoelectric power generator attached. 


thermometer (DE305) and compared with the measurement from 
our National Instrument PXIE4353 thermocouple module. 

The thermoelectric module as shown in Fig. 2(a) and (b) is a 
40(W)x40(L)x3.8(T) mm square diaphragm (GM250-127-14-16). 
Its technical data is shown in the top row of Table 1. The piezo¬ 
electric module is a circular diaphragm (Piezo Systems Inc. T216- 
A4NO-573X) made of Lead-Zirconate-Titanate, as shown in 
Fig. 2(c) and (d). It has thickness of 0.41 mm at resonant frequency 
of 290 Hz. It’s technical data is shown in the bottom row of Table 1. 

3. Experimental results 

3.1. Thermodynamic measurements 

To measure the temperature distribution in the bifurcating 
branches, an infrared thermal imaging camera is used. It captures 
the quartz tube surface temperature spectrum. However, it can 
characterize the ‘hot’ status of the oscillating flow inside each 
branch. The camera is set to operating range of 0 to 500 °C. Fig. 3(a) 
shows the measured temperature contour of the present setup in 
comparison with the temperature measurement of a conventional 
Rijke-type thermoacoustic system as shown in Fig. 3(b). It can be 
seen that the present system is associated with two different 
temperature ‘exhaust’ flows. The surface temperature of the upper 
short branch is around 26 °C, indicating that the air flow inside the 
branch is very close to room temperature, while the other upper 
branch is associated with ‘hotter’ air flow. However, the conven¬ 
tional Rijke-type system involves with only ‘hot’ exhaust oscillating 
flow, which limits the application of a piezoelectric generator due 
to the high temperature. Thus our present thermoacoustic system 
is unique in terms of producing dual-temperature oscillations, and 
providing a platform to apply energy harvesters with both ther¬ 
moelectric and piezoelectric energy conversion processes. 

To validate the temperature measurement, two arrays of K-type 
thermocouples are implemented on the bifurcating tube: one array 
with 4 thermocouples connected to the longer upper branch of 
length 0.9 m, the other with 3 thermocouples attached to the 
shorter branch of length 0.5 m. The measured temperatures are 
summarized in Fig. 4. It can be seen that one branch is with a ‘hot’ 
flow, while the other one is at ambient temperature. These mea¬ 
surements are consistent with the infrared camera measurement. 
The measurement difference between the infrared camera and the 
thermocouples is approximately ±2 °C. 

The presence of the thermoacoustic oscillations can be 
confirmed by measuring the pressure fluctuations along the 
bifurcating branches, since they cause the air pressure fluctuating. 
For this, two arrays of microphones (B&K Type 4954) are imple¬ 
mented and placed side by side to the thermocouples. The mea¬ 
surements are shown in Fig. 5(a) and (b). It can be seen that there 
are intensive pressure oscillations in both bifurcating branches, and 
the maximum sound pressure level is about 137 dB. The dominant 
mode of the thermoacoustic oscillation is at around 180 Hz. And it 
has several harmonics (see Fig. 5(a) and (b)), which indicates the 
nonlinearity of the system. This frequency corresponds to the 
fundamental mode with a wavelength twice of the total length of 
the bottom stem and the upper longer branch, while a conventional 
Rijke-type system with total length of 0.8 m is associated with 
around 230 Hz oscillations, as shown in Fig. 5(c). The ‘oscillating’ 
nature of thermoacoustics in the bifurcating branches is further 
confirmed via direct measurement by using a HWA (hot wire an¬ 
emometers) near the exit (about 0.5 cm) of the ‘cold’ branch, as 
shown in Fig. 5(d). This confirms that the present setup is a reliable 
system to produce dual-temperature oscillations. 

It is worth noting that the angle a between the two bifurcating 
branches plays an important role in producing dual-temperature 
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thermoelectri module 



Fig. 2. Photo of thermo- and piezo-electric power generators and their implementation. 


oscillations. To confirm it, temperature contour measurement of 
another bifurcating tube with a = 90° is conducted, as shown in 
Fig. 6(a). It can be seen that both ‘hot’ oscillating flows are pro¬ 
duced. The oscillatory nature of the ‘exhaust’ flow can be observed 
in the measurement of the pressure fluctuation phase diagram 
along the left and right branches, as shown in Fig. 6(b). The ‘hot’ 
oscillations produced are due to the buoyancy generated, which is 
related to the height of the bifurcating branch. If the buoyancy force 
is characterized by reduced gravity g [21], then it can be shown as 

g'=gn ^ (l) 

where H denotes the height of the upper branch open end with 
respect to the bottom open end. When a = 90°, g' L = g^. This gives 
rise to both bifurcating branches are involved with ‘hot’ oscillating 
flows. In order to study the effect of the angle a at different values 
on producing dual-temperature oscillations, further investigations 
are needed to determine the optimum geometric configuration. 

3.2. Electrical power generation measurement 

The oscillations consisting of both ‘hot’ and ‘cold’ pressure 
fluctuations might have potential application in industry, such as 
energy efficient drying for small-size particles such as soya beans or 


Table 1 

Parameters of thermo- and piezo-electric power generators. 


Couples 

> 

U 

fmax (A) 

P max (W) T max (C) R in (O) 

127 

9.4 

1.09 

5.10 

250 

7.9 

Weight 

Stiffness 

Capacitance 

Blocked 

Diameter 

Free deflection 

(g) 

(N/m) 

(nF) 

force (N) 

(mm) 

(pm) 

9.8 

5xl0 3 

107 

±2.4 

63.5 

±476 


coffee beans [6]. However, in the present work, we demonstrate 
another potential of using such oscillations for energy harvesting. 
This is achieved by implementing a thermoelectric generator in the 
‘hot’ bifurcating branch and a piezoelectric generator in the branch 
with ‘cold’ oscillations, as shown in Fig. 1. Thus the thermal energy 
harvested is maximized by simultaneously using two different 
energy conversion processes, which is different from the conven¬ 
tional standing-wave thermo-acoustic-piezo system 18]. 

Fig. 7 shows the electric voltage output and the power gener¬ 
ated from both the TE (thermoelectric) and PE (piezoelectric) 
modules. It can be seen that the total power output from PE and TE 
is about 5.71 mW, of which about 5.5 mW is produced by the TE. 
The PE generates approximately 0.21 mW electricity, which is much 
lower than TE. However, when it is compared with the electrical 
power generated from a similar PE implemented on a conduction- 
driven thermoacoustic system [18 , our PE generates 61% more 
power, as shown in Fig. 7(d). The PE performance is related to the 
oscillation frequency. In order to maximize its electricity output, it 
is tuned to ensure that its resonance frequency is closed to the 
oscillation frequency 180 Hz by adding an aluminum disk-plate 
weighting 2.85 g at the center of the diaphragm. Similar tuning 
was conducted by Smoker et al. [18]. Note that the electrical power 
output E l is calculated by: 

El = V( ms /R L = V 2 /2 R l (2) 

where l/ rms = V/V2 denotes the root mean square of the output 
voltage V. 

The energy conversion efficiency from thermal to electric power 
is also measured. Fig. 8 shows the overall efficiency of converting 
input heat power into output electrical power via the thermo¬ 
electric rie and piezoelectric modules ri a . It can be seen that the 
maximum ri e is about 0.0013%, when the loading resistance is 
approximately 7 Q. The overall energy conversion efficiency of the 
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(a) 



Temperature 
(500.0) °c 



( 0 . 0 ) 



Fig. 3. Temperature contour measurement by using infrared thermal imaging camera (a) dual-temperature thermoacoustic system (b) conventional Rijke-type system. 


TE module p e is much larger than that of PE p a . However, when the 
present system is compared with the thermo-acoustic-piezo sys¬ 
tem as reported by Smoker et al. [18 , the maximum p a is found to 
increased by 60%, as showed in Fig. 8(b). 

In order to gain insight on the heat-to-sound-to-electricity 
conversion, the conversion efficiency ri c from acoustical to elec¬ 
trical energy is measured. ri c is defined as ri c = E L /pV, where p and 
V denote the pressure and volume velocity at the piezo-diaphragm 
position. It can be seen that the piezoelectric diaphragm gives rise 
to the maximum energy conversion efficiency of approximately 
22% when the load resistance is 5070 Q. It is much greater than the 
maximum acoustical energy conversion efficiency of 9.5% obtained 
by Smoker et al. 18]. Furthermore, our system uses two different 
energy conversion processes. And it involves no heat exchangers 
and stack, which are critical components in conventional thermo- 
acoustic-piezo system 18]. Note that further investigations are 
needed to minimize the size of such dual-temperature thermoa¬ 
coustic system and to maximize its performance so that it can be 
practically applied. 

4. Theory of heat-driven Rijke-Zhao oscillations 

4 A. Thermodynamic analysis 


parameters, which are set to 0.1 and 1.65 respectively [23 . If ztf 
denotes the energy release per gram of the gas mixtures by the 
complete reaction, the rate at which heat is generated in the gas per 
gram is given by 

Q(T, 0 = ~ [C 3 Hs] t _ 0 d[( ^ sl = Q(T) + Q!(T, t) (6) 

where the instantaneous heat release Q(T,t) consists of a mean 
value Q.(I) and a perturbation part O(T,t). 

If we denote the internal energy of the gas mixture by E, then 


CvP _ 1 p 

gfcp y - 1 p 



where c v is the constant-volume specific heat and M = ^ u /^w- 
The first law of thermodynamics considering energy conservation 
can be described by 


dE dQJJ\t) d/(l/p) 

dt at ^ dt 



The convection flow in the present system has negligible Mach 
number M a [24,25 , so that the flow properties in a region of the gas 


The pressure of the gas mixtures used in the present setup is al¬ 
ways low enough to justify the assumption of the perfect gas laws. 
Thus ifp denotes the instantaneous pressure, T the temperature and p 
the density, 


P 



JZ W 




where c is sound speed, y = 1.4 is the ratio of specific heats, &7 U is 
the universal gas constant (=8315 J/kmol-I<), is the molar mass 
of the gas mixture. Since the fuel burned is LPG (liquefied 
petroleum gas) with main component of C 3 H 8 22], the complete 
chemical reaction in single step Arrhenius kinetics is of the form 

C 3 H 8 + 50 2 ^3C0 2 + 4H 2 0 (4) 

Then the rate of the reaction is given as 

= -^ f e-^/^ T [C 3 H 8 ] b [0 2 ] c (5) 

where is the Arrhenius constant describing the collision fre¬ 
quency of molecules. is the activation energy. [•] denotes the 
concentrations of reactant species, b and c are the reaction rate 
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Fig. 4. Temperature measurements by using two arrays of K-type thermocouples along 
the bifurcating branches. 
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Frequency (Hz) 



Fig. 5. Measured thermoacoustic oscillations: (a) frequency spectra along the upper long branch, (b) frequency spectra along the upper short branch, (c) frequency spectra along the 
conventional Rijke-type system with total length of 0.8 m and (d) velocity fluctuation near the upper short branch open end. 


mixture where the mean pressure, temperature, density and ve¬ 
locity are p, T, p and u respectively can be expanded as 

rp 

|=l+M a <T 1 + C»(Mf<T 2 ) + --- ==l+M a T 1 + o(Mft 2 ) + - (9) 

r i 


^=1+M a X\+o(Mlx2)+--- ==l+M a /r 1 +e>(M al u 2 )+--- (10) 

Substituting the linearized form of Eqs. (9) and (10) by 
neglecting the second and higher order terms into Eq. (8) leads to 


d 

dt*“l v 


ln{(l+M a T 1 )(l+M a x 1 ) 1 r 


T - 1 


aQ(T, t) 


(i+M a T!)^r at 


(ii) 

Application of Taylor series expansion on Eq. (11) and 
neglecting second order terms gives rise to the linearized form of 
Eq. (11) as follow 


y — 1 dQ(T, t) a ^ , 


( 12 ) 


The motion of the perfect gas can be described by the following 
equation 


where the vector V denotes the gas velocity. And it can be shown 
that dV/dt is given as 


dv av av _/1._ 2 


dt 


at 


+ (V-V)V 


at + V V2 V ' -Vx(VxV) 


(14) 


Since the gas flow under consideration is irrotational i.e. 
curlV = 0, a velocity potential 4> exists [26] as follow 


V = -V<p 

Substituting Eqs. (3), (14) and (15) into Eq. (13) gives 


M 


vl 5 v2 


_ mr 

= mvr + —vp 


dt \Z ) p 

The equation of mass conservation can be shown 


V'0>V) + ! 


0 


dt 


Substituting Eqs. (9), (10) and (15) into Eq. (17) yields 


{ln(l +M a X!)} = W 2 <p 


(15) 


(16) 


(17) 


(18) 


dV 

dt 


1 

+ -Vp 
P 


0 


(13) 


1 When Q(T,t) = 0 (i.e. no heat source), Eq. (11) can be simplified to 
(1 + M a T -[) = (1 + MaXi ) T_1 > which describes the general adiabatic law in terms of 
M a X! and M a 


Use of Taylor series expansion and neglecting second and high 
order terms leads to 


a ( MaX b = v 2 <f> 

dt 


( 19 ) 
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Fig. 6. (a) Temperature contour measurement of a bifurcating tube with 2a; = 180° (b) 
phase diagram between p(t ) and p(t) along the left and right branches. 


Similarly, we can linearize Eq. (16) by neglecting second order 
terms of M a T\ and M a Xi and such products as M a T! VA/f a Xi, and 
integrating it, then we obtain 


However, when there is thermal energy input (i.e. 
50.(7, t)/0t=£ 0), Eq. (21) can be reformulated as 


d 2 (p 

0? 


c 2 v 2 0 = (Y-l)^Q(r) = (Y-t)£ Q(T,t)+Q'(T,t) 


d 


d r- 


dt 


at 


(23) 


Eq. (23) describes the heat-driven acoustic wave equation. 
Further analysis on the right-hand-side of Eq. (23) shows that the 
unsteady heat addition Q(7,t) plays a dominant role in producing 
acoustic waves, while the mean part Q(7) does not make contri¬ 
bution. Q(7) is convected away by using heat exchangers as typi¬ 
cally found in the conventional standing-wave system 18]. 
However, if Q(7) is harvested by using thermo-electric power 
generation module, then the thermal energy harvesting is maxi¬ 
mized. This is achieved in our experiment. 

If the acoustic damping effect such as the boundary losses, 
vorticity and flame-holders induced drag is involved, Eq. (23) can 
be shown as, 



(r-D^Q(T,t) 


(24) 


where £ denotes the overall damping coefficient. Substituting Eq. 

(6) into (24) leads to 



C 2 V 2 </> 


(1 -Y)^9 2 [C 3 H 8 ] 
[C 3 H 8 ] t=0 at 2 


(25) 


Since thermodynamic parameters such as pressure or density 
can be expressed in terms of velocity potential 4>, Eq. (25) describes 
a general nonlinear heat-driven acoustic oscillation formulation by 
burning C 3 H 8 . The heat-driven acoustic oscillation equation can 
also be applied to other fuels. 


4.2. Thermoacoustic model 


In order to illustrate the generation mechanism of acoustic os¬ 
cillations by unsteady heat input and to qualitatively compare with 
experimental measurements, numerical simulations are conduct¬ 
ed. A bifurcating tube with a localized heat source confined is 
simulated, corresponding to the experimental setup as shown in 
Fig. 9. The momentum and energy equations across the heat source, 
when linearized with respect to the perturbation quantities, can be 
written as 


1 dp' _ 0 u' 

?dF + p aT 


(7-1) 

7=2 



(26) 


^=«r(M a t 1 +M a x 1 ) (20) 

Differentiating Eq. (20) with respect to t and replacing 

0(M a i 1 )/at and 0(M a Xi)/6t with Eqs. (12) and (19) respectively 
gives 

= (Y — l)^Q(7\t) (21) 

where the speed of sound c 2 = yMT. When there is no thermal 
energy transferred into the system (i.e. Q(7,t) = 0), Eq. (21) reduces 
to the well-known acoustic wave equation 

^-c2 V 2 0 = 0 (22) 


2 Note that any additional term which depends on the time and not on the space 
coordinates is neglected from 4>. 


du' 



dp' 

dX 


0 


(27) 


By differentiating Eq. (26) with respect to time t, Eq. (27) with 
respect to x and subtracting, an inhomogeneous heat-driven wave 
equation is obtained namely, 


1 d 2 p' £ dp' d 2 p' 

?"dt r + ?"aF - a* 2- 


y - 1 0 

C 2 at 


(Q'(t)<5(x-x f )) 


(28) 


where <5(x) is dirac delta function describing the localized heat 
source. 

Expanding the pressure perturbation as Galerkin series [4,27] 
gives, 

N i 

P'M) = E —4(W m W (29) 

1 rim 

m = 1 111 111 

where K m is a coefficient. Its physical meaning and derivation are 
shown later. N denotes the mode number. Overdot denotes the time 
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(a) 



(c) 



Fig. 7. Measured voltage and power output from the thermo 


(b) 



(d) 



and piezo-electric generators with varying resistance in Q. 


derivative. The function ^ m (x) are the eigen-solutions of the ho¬ 
mogeneous wave equation, d 2 p'/dt 2 - c 2 d 2 p'/dx 2 = 0. 

The presence of the bifurcating branches at x d , and the heat 
source resulting in the mean temperature undergoing a jump from 
the region 0 < x < x f to x f < x < x d give rise to 4 different flow 
regions. Substituting Eq. (29) into the homogeneous wave equation 
gives 



0<x<x f 


h 


xsin ^ ,Jm( ^ 


+cos( ^ m( j 


Xf <x<x d 




h 


sin 


f M m (x d -P 


sin 


d3 , 




/ 6Jm(X~X d3 ) \ 

y c 3 cosa; ) 


X d <x<x d3 



d 2 jm 

dx 2 


0 


(30) 


h 


sin 


( (x d -; 


sin 


•d4. 




/ 6Jm(X-X d4 ) \ 
\ C4COSQ; ) 


x d <x<x d4 


where q denotes the sound speed at ith region. It is assumed that 
there is no mean temperature jump from the region x f < x < x d to 
x d < x < x d4 . Thus c 4 = c 2 and u 4 = u 2 . Mass and pressure conti¬ 
nuity hold at x = x d as P 2 U 2 S 2 + P 3 U 3 S 3 = p 4 u 4 S 4 and 
P 2 (x d ,t) = P 3 (x d ,t) = p 4 (x d ,t). After linearizing them and neglecting 
the higher order terms, two jump conditions at x d are given as 


where h and K are given as 

c 2 S 4 cot(( dm (x d - x d4 )/c 4 cosa) 
S 2 c 4 cosa 

c 2 S 3 cot(a) m (x d - x d 3 )/c 3 cosa) 
S 2 c 3 cosa 


(32) 


(33) 


= 0 and [^ m ]| g = 0 (31) 

Applying the pressure continuity and velocity jump across the 
heat source and the boundary conditions at x = 0, x = x d 3 and x = x 4 
[4,27], i.e. au'(x, t)/dx \ x=0 = 0, au'(x, t)/9x| x=Xd3 = 0 and 
0u'(x, t)/ax| x=Xd4 = 0 leads to 


dx 


h 


sinfwXf/c! 


xsin(w m (x f - x d ) j c 2 ) + cos(w m (x f - x d ) j c 2 ) 


(34) 


where S denotes the cross-sectional area of each branch. 
0 < a < 90° is half of the angle between the two bifurcating 
branches. 

i// m (x) satisfies the condition of orthogonality, i.e. 
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Fig. 8. Measured overall energy conversion efficiency from the present thermoacoustic 
system with the thermoelectric generator and tuned piezoelectric diaphragm. Note 
that the ‘tuned’ and ‘untuned’ denote the aluminum disk ‘attached’ and ‘unattached’ 
respectively. 


*d *d4 
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0 *d 

Further analysis illustrates that the eigenmode frequency oj m is 
satisfied with 


, . xs in( ^fa) )+cos 




sin 


, (x f —x d ) 


C 2 


(36) 


Substitution for p (x,t) from Eq. (29) into Eq. (27) lead to 


£:(^ + u 2 m V m y m (x) = (T-l)^«5(x-x f ) (37) 

After multiplication by \// m (x) and integration with respect to x, 
the orthogonality condition of Eq. (35) shows that Eq. (37) can be 
simplified to 
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Fig. 9. Schematic of a dual-temperature thermoacoustic system with the thermo¬ 
electric and piezoelectric power generator attached. 
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Substituting Eq. (32) in Eq. (29) yields 


(39) 


Vm (0 

Km OJm 




iMl 
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p'(z,t) = < 


xsin( 6Jm( ^ Xd) ) +cos( ^ m( ^ * d) 


0<x<x f 


x f <x<x d 


( N 

E 

m = l 
N 

E 

m =1 

E i£hsin^^)/sin^§^) x d <x<x d3 

m = 1 x / / \ / 

N 


E ^sin(^|^)/sin(^^) x d <x<x d4 

(40) 

Substituting the Galerkin series expanded p (x,t) into Eq. (28), 
multiplying with i// m (x) on both sides and integrating from x = 0 to 
x = x d 4 yields 


k/n d 7J ~ K m d 7) m 
dt 2 dt 


+ K m (x) m 7) m 


-(y-l)Q'(t)sin(^) 

(41) 


The system Eq. (41) can solved if Q.(t) is known. In our experi¬ 
ment, the heat source is a premixed laminar flame. However, for 
simplicity, the localized heat source is numerically modeled as an 
electrical heater consisting of heated cylindrical hot wire in a 
perpendicular flow. It unsteady heat release rate is described by 
using a modified form of King’s Law [29] as: 



where n m is given as 


where L w is the length of the wire, T w is the wire temperature, T is 
the temperature of the surrounding air, x is the heat conductivity of 






































































































D. Zhao et al. / Energy 65 (2014) 517—526 


525 


the air, c v is the specific heat of the air per unit mass at constant 
volume, p is the mean density of the surrounding air and d is the 
diameter of the wire and Si is the cross-sectional area of the tube. u' f 
is the oncoming flow velocity and it can be estimated by using Eqs. 

(27) and (40) as 



N 


E 

m = 1 


Vm 

Pl Ci K m 


cos 



(43) 


Now the time evolution of the heat-driven acoustic perturbations 
from specified initial conditions can be determined by numerically 
solving Eqs. (41) and (42) using a 4th Runge-Kutta method. In our 
calculations, the duct radius is 0.025 m, the angle between the two 
bifurcating branches with lengths ofxd 3 = 0.783 and Xd 4 = 1.169 m is 
a = 15°, the length of the mother tube is Xd = 0.3 m, and the heat 
source is placed at Xf = 0.13 m away from the bottom open end. The 
mean heat release rate Q.(T) is set to 2.43 x 10 5 mW/m 2 . 


43. Rayleigh index and energy redistribution 


Rayleigh index ®. Following the work 28 , we define a normalized 
Rayleigh index ©(£) as given as 



It can be seen that when initial flow disturbances grows into a 
limit cycle, unsteady heat release is in phase with the pressure 
fluctuation, i.e. @(w l5 t) > 0. This phase information is consistent 
with Rayleigh’s criterion 7], verifies our numerical model and 
explains why self-sustained thermoacoustic oscillations occur in 
the bifurcating tube. 

To gain insight on the energy variation of the thermoacoustic 
modes, the energy redistribution between various modes are esti¬ 
mated. Fig. 11 shows the energy variation of the first 4 modes with 
time. The acoustical energy (E(t)) at is defined as 


Fig. 10(a) shows that the self-sustained thermoacoustic oscilla¬ 
tion in terms of pressure is generated by introducing an initial 
impulse perturbation. The initial disturbance is quickly grown into 
limit cycle, as shown in Fig. 10(a). The resulting oscillations are at 
around 180 Hz with a sound pressure level 138 dB. Harmonics as an 
indicator of nonlinearity are also observed, as shown in Fig. 10(b). 
Comparing these numerical results with our experimental mea¬ 
surements, good agreement is obtained. Fig. 10(c) illustrates the 
variation of Rayleigh index @ with time, which is generally used as 
a critical indicator of the thermoacoustic system stability [7,28]. It 
describes the energy exchange rate between unsteady heat release 
and the pressure oscillations. If ® > 0, the acoustical energy in the 
present system is increased and tends to drive the unstable process 
into a limit cycle. However, a negative @ indicates a ‘destruction’ 
interaction between unsteady heat release and sound, and so the 
flow oscillations decay. Thus it would be interesting to check the 
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Fig. 10. Heat-driven acoustic oscillation in a bifurcating tube: (a) a limit cycle, (b) 
frequency spectra, (c) Rayleigh index 0 1 (w 1 , t) variation with time. 
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(45) 


The initial triggering disturbances with a small amplitude, such 
as impulse noise is set to 6 ) 4 . However, the fundamental mode with 
the lowest frequency w 1 is excited and grown into a limit cycle. As 
the initial disturbances are decaying, thermoacoustic modes at 6)2 
and 6)3 also gain energy. The energy redistribution might result 
from the direct interaction between the two eigenmodes interact 
directly, or the interaction of various modes with the base flow. 





Fig. 11. Sound energy re-distribution: (a) variation of log 10 ((,£i (w \, t),}) with time (b) 
variation of log 10 ((,£ 2 (^ 2 , t)>)) with time (c) variation of log 10 ((,£ 3 ( 61 ) 3 , t),)) with time 
(d) variation of log 10 ((,£ 4 (w 4 , t),)) with time. 















































































526 


D. Zhao et al. / Energy 65 (2014) 517—526 


5. Conclusions 

Thermodynamic measurement and analysis of dual¬ 
temperature thermoacoustic oscillations for energy harvesting 
application are conducted. A bifurcating tube is designed and tested 
to generate both ‘hot’ and ‘cold’ thermoacoustic waves. The acoustic 
fields in the bifurcating tube are monitored by implementing mi¬ 
crophones, thermocouples, a hot wire anemometer and an infrared 
camera. The maximum sound pressure level is found at around 
138 dB. Energy harvesting is achieved by using two different energy 
conversion processes via implementing both piezo- and thermo¬ 
electric power generation modules. This is different from the con¬ 
ventional standing-wave-piezo system [18 , of which heat ex¬ 
changers are involved. And only acoustical energy is harnessed, 
while the ‘waste heat Q.’ is dumped. It is found in our experiment 
that the total electric power output is 5.71 mW, of which approx¬ 
imately 5.5 mW is produced by the thermo-electric module. The 
power output is much larger than that harnessed from the 
standing-wave-piezo system. 

To gain insight on the generation mechanism of the dual¬ 
temperature oscillations, a ‘T-shaped’ tube is tested. It is found 
that the angle between the bifurcating branches plays an important 
role. The angle effect is then considered in developing a nonlinear 
model to simulate our experiment. Thermodynamic laws are first 
used to derive a general formulation to illustrate how acoustic os¬ 
cillations are produced by inputting thermal energy via burning 
C 3 H 8 . The formulation is shown to be applicable to other type of 
fuels. It is also shown that unsteady heat release is an efficient 
sound source. And nonlinearity characterized by harmonics is 
associated with the thermoacoustic system. In comparison with our 
experimental measurement, qualitatively good agreement is ob¬ 
tained in terms of the oscillation frequency and sound pressure 
level. Rayleigh index ® is then examined to gain insight on the 
interaction between unsteady heat release and oncoming flow 
disturbances. It is found that when @ > 0, sound energy converted 
from thermal energy keeps increasing and a limit cycle is eventu¬ 
ally generated. Finally, the energy redistribution between various 
eigenmodes is examined. Sound energies of higher frequency 
modes are found to be transferred to the lower frequency ones. And 
lower frequency modes are easier to trigger. 
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